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I. INTRODUCTION 



The (/-state Potts antiferromagnet (AF) |],|2j exhibits nonzero ground state entropy, So > 
(without frustration) for sufficiently large q on a given graph or lattice. This is equivalent 
to a ground state degeneracy per site W > 1, since So = kshiW. Such nonzero ground 
state entropy is important as an exception to the third law of thermodynamics ||. The zero- 
temperature partition function of the above-mentioned g-state Potts antiferromagnet on a 
graph G satisfies Z(G,q,T = 0)paf = P(G,q), where P(G,q) is the chromatic polynomial 
expressing the number of ways of coloring the vertices of the graph G with q colors such that 
no two adjacent vertices have the same color ||- 0. Thus 

W({G},q) = lim P{G,q) 1/n (1.1) 

ra^oo 

where n = v(G) is the number of vertices of G and {G} = lim^oo G. The minimum number 
of colors needed for this coloring of G is called its chromatic number, x{G)- At certain 
special points q s (typically q s = 0, 1, .., x(G)), one has the noncommutativity of limits 

lim lim P{G,q) 1/n ^ lim lim P{G,q) 1/n (1.2) 

q— >q s n— +oo n— >oo q—>q s 

and hence it is necessary to specify the order of the limits in the definition of W({G},q s ) 
0. As in 0, we shall use the first order of limits here; this has the advantage of removing 
certain isolated discontinuities in W. In addition to Refs. 0- 0, some other previous related 
works include Refs. ||- p2| . Since P(G,q) is a polynomial, one can generalize q from Z + 
to R and indeed to C. W^({G}, q) is a real analytic function for real q down to a minimum 
value, q c ({G}) |7|,p^|. For a given {G}, we denote the continuous locus of non-analyticities of 



W as B. This locus B forms as the accumulation set of the zeros of P(G, q) (chromatic zeros 
of G) as n — > oo ||,|l(],|r|,0] and satisfies B(q) = B(q*). A fundamental question concerning 
the Potts antiferromagnet is the form of this locus for a given graph family or lattice and, in 
particular, the maximal region in the complex q plane to which one can analytically continue 
the function W({G},q) from physical values where there is nonzero ground state entropy, 
i.e., W > 1. We denote this region as R±. Further, we denote as q c the maximal point where 
B intersects the real axis, which can occur via B crossing this axis or via a line segment of 
B lying along this axis. 

In the present work we present exact calculations of P(G, q) and ^({G}, q) for families of 
cyclic graphs consisting of linked polygons, where the polygons may be adjacent or separated 
by a certain number of bonds. By the term cyclic we mean that these families of graphs 
contain a global circuit, defined as a route along the graph which has the topology of S l and 
a length £ g . c . that goes to infinity as n — > oo. The results provide new insight into W and B 
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in an exactly solvable context where one can study their dependence on several parameters 
characterizing the families of polygon chain graphs. The present work extends our recent 
studies of cyclic families of graphs [f7tp2"l|. 

From our previous exact calculations of P(G, q) and W({(j}, q) on a number of families 
of graphs we have inferred several general results on B: (i) for a graph G with well-defined 
lattice directions, a sufficient condition for B to separate the q plane into different regions is 
that G contains at least one global circuit |I9|P1, and (ii) the locus B for such a graph does 



not contain any endpoint singularities. Two other general features are that for graphs that 
(a) contain global circuits, (b) cannot be written in the form G = K p + H f\ and (c) have 
compact B, we have observed that B (iii) passes through q = and (iv) crosses the positive 
real axis, thereby always defining a q c . The families of polygon chain graphs studied in this 
paper do contain global circuits, and the exact results presented here are in accord with the 
above general properties. 

The chromatic polynomial P(G, q) has a general decomposition as 

N a 

P(G, q ) = c ( q )+J2cMM<l)) tjn (1.3) 

where the aj(q) and Cj^{q) are independent of n, while Co(q) may contain n-dependent terms, 
such as (— l) n , but does not grow with n like {const.) 11 with \const.\ > 1. The expression 
Co may be absent. A term ai(q) is "leading" if it dominates the n — > oo limit of P(G,q). 
The locus B occurs where there is an abrupt nonanalytic change in W as the leading terms 
an changes; thus this locus B is the solution to the equation of degeneracy of magnitudes of 
leading terms. Hence, W is finite and continuous, although nonanalytic, across B. 
It is convenient to define the following polynomial: 

B *(o> = §pff = «'- 2 £(-^ = D-D'f ~ < L4 > 



1 Some families of graphs that do not have regular lattice directions have noncompact loci B that 
separate the q plane into different regions [16,18,2 1 1 . 



2 The complete graph on p vertices, denoted K p: is the graph in which every vertex is adjacent to 
every other vertex. The "join" of graphs G\ and G2, denoted G\ + G2, is defined by adding bonds 
linking each vertex of G\ to each vertex in G2. 
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where 



a = q-l (1.5) 

and P(Ck, q) is the chromatic polynomial for the circuit (cyclic) graph Ck with k vertices, 

P(C k ,q) = a k + (-l) k a (1.6) 

This paper is organized as follows. In Sections II and III we present our results for P, W, 
and B for the open and cyclic polygon chain graphs. Some concluding remarks are included 
in Section IV, and some further results are included in the Appendix. 



II. OPEN POLYGON CHAIN GRAPHS 

Before giving our results for the cyclic polygon chain graphs, we discuss the simpler case 
of families of open polygon chain graphs comprised of m subunits, each subunit consisting 
of a p-sided polygon with one of its vertices attached to a line segment of length e g edges 
(bonds). Thus, the members of each successive pair of polygons are separated from each 
other by a distance (gap) of e g bonds along these line segments, with e g = representing the 
case of contiguous polygons. Since each polygon is connected to the rest of the chain at two 
vertices (taken to be at the same relative positions on the polygons in all cases), the family of 
graphs depends on two additional parameters, namely the number of edges of the polygons 
between these two connection vertices, moving in opposite directions along the polygon, e\ 
and e 2 . For e\ ^ e 2 we denote the smaller (s) and larger (£) edge lengths as 

e s = min(ei, e 2 ) , = max(ei, e 2 ) (2.1) 

Clearly, 

e 1 + e 2 =p (2.2) 

We denote this family of open (o) polygon chain graphs as G ei ^ 2>eg>m - . Illustrations of these 
graphs are given in Fig. ([I]). Since it is immaterial whether one labels the upper and lower 
routes between the left and right-hand points on each polygon where it connects to the rest 
of the chain as t\ and e 2 or in the opposite order, as e 2 and e\, 

G^ei,e 2 ,e g ,m;o = C e2jeij e 9jm;0 (2-3) 
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FIG. 1. Illustrations of cyclic and open polygon chain graphs G ei , e2 , e m and G eij e 2j e 3 ,m;o with 
(d, e2, e g , m) = (a) (2,2,0,6), (b) (2,2,1,4), (c) (2,3,2,3). For the cyclic (open) graphs, the rightmost vertex 
on each graph is identified with (is distinct from) the leftmost vertex at the same level, respectively. 

The number of vertices is 

u ( G ei ) e a ,e fl) m;o) = (p + e g - l)m + 1 (2.4) 

and the number of edges or bonds is 

e(Gei,C2 ) e fl> m;o) = (P + Cg)™ (2.5) 
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The chromatic number x — 2 if p is even and x = 3 if p is odd. The chromatic polynomial 
is easily calculated to be 



P(G ei , e2 , eg , m . ,q) = q(a l ) m (2.6) 

where 

a x = {q- l) e <> +1 D p (2.7) 

Evidently, the chromatic polynomial has the form of ( |1.3|) with cq = 0, N a = 1, c\ = q{ai)^ tl 
and ti = t, where 

t = (2.8) 

P + e g - 1 

Note that this chromatic polynomial does not depend on ei or e 2 individually, but instead 
only on their sum p, i.e., the number of sides of each polygon^. 

Let us next consider the limit of an infinitely long chain, i.e., m — > oo. Formally, we 
denote this limit as 

L m : m — > oo with e\ , e-i , e g fixed (2.9) 

We find 

W([ lim G ei , e2<eg , m . ],q) = [(q-l) e " +1 D p )^ (2.10) 

The nonanalyticities of this W function consist of discrete branch point singularities at its 
zeros, including the point q = 1 and the p — 2 zeros of D p ; there is no continuous locus of 
nonanalyticities of W\ i.e., B — 0. If p is odd, then D p contains a factor (q — 2), so that one 
of the zeros of W is at q = 2. 

Given the relation (iO), there are also other ways to take the number of vertices to 



infinity: 

L eg '■ e 9 — * oo with e\ ,e 2 ,m fixed (2-H) 



3 One can consider a generalization of this family in which the quantities e\ and e 2 can be different 
for each polygon. The chromatic polynomial and W function for this much larger family of graphs 



are again given by (0) with <SM and by pT0|), respectively. 
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L ei : e\ — > oo with e 2 , e g , m fixed 



(2.12) 



-^e 2 : e 2 — ► oo with e\ ,e g ,m fixed (2-13) 

and 

L p : p — >■ oo with e2 — e± ,e g ,m fixed (2-14) 
For the L e limit, we find 

W([ lim G ei , e2 ,e 9 , m;0 ],g) = g- 1 (2.15) 

e g — >oo 

This is analytic for all q so that, in particular, B = 0. The W function in eq. (|2.15|) is the same 
as the W function for linear and, more generally tree, graphs. This is understandable, since 
as e s — ■> oo, the vertices on the linear connecting segments occupy a fraction approaching 
unity of the total number of vertices of the graph. 

Since the chromatic polynomial (|2.6|) only depends on e\ and e 2 through the quantity p, 
it follows that for this family of open polygon chain graphs, the limits L ei , L e2 , and L p all 
yield the same result. Further, from eq. ( [L.4|) , it follows that in these limits, W is the same 
as that for a circuit graph 0, viz., 

W([\imG eue2 , egim , },q) = q-l for |g-l|>l (2.16) 

p— >oo 

andQ 

|^([limG ei ,e 2 , e9 , m; o],g)| = l for |g-l|<l (2.17) 

For this limit, B consists of the circle \q — 1| = 1. Again, these results are easily under- 
stood, since as p — » oo, the vertices on the polygonal circuit subgraphs occupy a fraction 
approaching unity of the total number of vertices. 



4 For real < q < q c ({G}), as well as other regions of the q plane that cannot be reached by 
analytic continuation from the real interval q > q c ({G}), one can only determine the magnitude 
|W({G},<7)| unambiguously jjj. 
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III. CYCLIC POLYGON CHAIN GRAPHS 



We proceed to our main subject, families of graphs consisting of m p-sided polygons 
linked together to form closed circuits. The other parameters, e g , e\, and e 2 , are the same as 
for the open polygonal chains discussed above. Thus, a given cyclic polygonal chain graph is 
specified by the four parameters (ei, e 2 , e g , m), so we denote it as Gr ei ,e 2 ,e 9 ,m> where the cyclic 
property is implicitly understood. For the same reason as given above for the open chain, 
we have 

The number of vertices in the graph G eii e 2 ,e 9i m is 

n = V (G ei ,e 2 ,e gi m) = (P + £g ~ I)™ (3.2) 

and the number of edges (bonds) is the same as for the open chain, e(G eiie2)egjm ) = (p + e g )m. 
If p is odd, then the chromatic number \ is 3. If p is even, then x = 2 or 3 depending on 
the values of e\, e 2 , e g , and m. In Fig. (H) we show illustrative examples of families of cyclic 
polygon chain graphs. From ( p.l|) , it follows that 



P(Ge 1 ,e2,e g ,m, <?) — P(Ge 2 ,e 1 ,e g ,m, <?) (3-3) 

Using the deletion-contraction theorem to get recursion relations which we then solve, 
we calculate the chromatic polynomial 

P(G ei , e2>eg>m , q) = (ai) m + (q- l)(a 2 ) m (3.4) 
where the term a± is the same as the term appearing in the chromatic polynomial for the 



open polygonal chain, ( |2.7p , and 

a 2 = (-ly+^q- 1 [g - 2 + (1 - g) ei + (1 - g) 62 ] (3.5) 



s=2 v 7 s=2 v 



The chromatic polynomial Q3.4|) has the form of ( |1.3|) with c = 0, A^ a = 2, c\ = 1, c 2 = (q— 1), 
and ti = t 2 = t, where t was given in eq. (|2~8|). Note that, owing to the second term, a 2 , the 
chromatic polynomial (|3.4| ) depends on the values of e\ and e 2 individually, in contrast to 
the result ( j2.6|) for the open polygon chain, which only depends on e± and e 2 through their 
sum, p. 
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If the smaller distance e s = 1, then the chromatic polynomial Q3.4|) reduces to the factor- 
ized form 



P(G es=1>ei , eg , m , q) = (D p ) m P(C nc ,q) = q(q-l)(D p ) m D nc (3.7) 

where 

n c = m(l + e g ) (3.8) 
Having calculated the chromatic polynomial (|3.4|), we can next obtain the W function and 



its continuous locus of nonanalyticities B for the various limits ( |2.9| ) and fl2.11| )-( p.l4j ). This 
locus is determined as the solution of the degeneracy equation of (nonzero) leading terms 

Kl = \a>2\ (3.9) 

It will be useful to reexpress this in the variable a = q — 1 of eq. ( |1.5| ). Multiplying both 
sides of ( J3.9|) by \a + 1|, we obtain 

\a e *{a p + (-l) p a)| = \a - 1 + (-a) £l + (-a) e2 | (3.10) 



(Since the value q = 0, i.e., a + 1 = 0, is already a solution of (|3.9| ), multiplying both sides 
of this equation by \a + 1| to get ( |3.10| ) does not introduce a distinct spurious zero.) 

We first discuss the limits that yield the simplest results. For the L Sg limit, B consists of 
the unit circle \q — 1| = 1 and divides the q plane into two regions, the exterior and interior 
of this circle, which regions are denoted as R\ and R2, respectively. If \q — 1| > 1, then ai 
is the dominant term and 

W{[\hn G eue2>eg)m ],q) = q-1 for |g-l|>l (3.11) 

e g — >oo 

If I q — 1 1 < 1 , then 02 is dominant and 

\W{[]xm G ei , e2 , e9 , m ],g)| = l for \q - 1| < 1 (3.12) 

That is, in the L g limit, the W function is the same as for the infinite- vertex limit of the 
circuit graph. This is expected, since as e g increases, more and more of the vertices of the 
graph lie on a circuit graph where, e.g., the route of the circuit follows either the e x or e 2 
bonds of a given polygon, then traverses the e g connecting bonds to the next polygon, and 
so forth around the cyclic chain. As a result of the symmetry ( |3.3|) , the L ei , L e2 , and L p 
limits are equivalent; in this case, we find that for \q — 1| > 1, the a\ term is dominant 
and W is again given by the right-hand side of ( |3.11|) while for \q — 1| < 1, neither term 
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dominates over the other, and \W\ = 1. The results for W are thus the same as for the L £g 
limit, and the reason is similar: as p — > oo, the vertices located on the polygons occupy a 
fraction approaching unity of all of the vertices. 

The most interesting and complicated results are for the L m limit. We have proved a 
number of general properties of B. We show calculations of B and comparison with chromatic 
zeros for long finite chains in Figs. 2-5 (and, in addition, Figs. 6-9 in the Appendix). Because 
of the factorization ( |3.7|) , the families with e s — 1 yield very simple results (see below) so 
that in the figures we concentrate on the cases with e s > 2: (ei, e 2 ) = (2, 2), (2,3), (2,4), and 
(3,3) and, for each pair (ei, e 2 ), the values < e g < 3. The general properties of B are 

• (Bl) B is compact. 

• (B2) B passes through q = 0. 

• (B3) B encloses regions in the q plane. 

• (B4) if e s = 1, then B is the unit circle \q — 1| = 1 independent of the values of tt 
and e g ; thus, in this case, q c = 2. The chromatic zeros of G ei>e2>e m lie exactly on this 
locus, independently of the values of ee, e g , and m. 

• (B5) if p is even, then q c = 2. 

• (B6) if p is odd and e s > 2, then q c < 2 and for fixed (e\, e<z), q c increases monotonically 
as e g increases, approaching 2 from below as e g — > oo 

To prove (Bl), we re-express eq. ( |3.9|) in terms of the variable 

y = - = 7 3.13 

a q — 1 



Then eq. ( |3.10|) becomes 



|1 + {-lfy p - l \ = \y e ^ es {y ei ~ x - y ee + (-1) ££ + (-l) es ^" es )| (3.14) 



where we have divided ( |3.9| ) by |a efl+p |, which factor is nonzero for q ^ 1 and thus, in 
particular, for the large-|g| (small \y\) region of interest here. The necessary and sufficient 
condition that B is noncompact, i.e. unbounded is that the equation (|3.14j) has a solution 



for y = 0. Clearly, it does not have such a solution, which proves property (Bl). 
For (B2), we use the property that 

D p (q = Q) = (-lY(p-l) (3.15) 
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so that |ai| = \p — 1| at q = 0. From the second equivalent form of a 2 in eq. (|3.6| ), it follows 
that |a 2 | = \p — 1| at q = also, so that the point q = is a solution to eq. (|3.9| ) and hence 
is on the locus B , which proves (B2). 

Property (B3) can be proved by explicit solution of (|3.9|) . One also observes from the 
results presented in Figs. 2-9 that in cases where B does not contain any multiple points^, the 
number of regions, N reg , and the number of connected components on B satisfy the relation 
N reg . = N comp , + 1. As was the case in our earlier work plf , one can get a general upper 



bound on the number N comp , using the Harnack theorem [^3] from algebraic geometry, but 
it is not very restrictive. To do this, we write the equation ( |3.9|) whose solution set is B 
out in terms of the real and imaginary components of q or, more conveniently, of a — q — 1. 
This yields a polynomial equation of general degree d = 2{p + e g ) in these components. 
In the cases without singular (multiple) points on B, the Harnack theorem then gives the 
upper bound N comp , < h + 1, where h, the genus of the algebraic curve comprising B, is 
h = (d — l)(d — 2)/2 p3| . Thus, for example, for the lowest case (ei,e 2 ,e g ) = (2,2, 1), the 
Harnack theorem gives the upper bound N comp , < 37. Clearly, this is a very weak upper 
bound compared with the actual result N comp , = 2. The upper bound becomes even higher 
for larger values of e g , and hence weaker, since the result remains iV comp = 2 for all families 
of the form (ei, e 2 , e g ) = (2, 2, e g ). Similar remarks apply for the Harnack bound applied to 
other cases (ei,e 2 ,e s ). 

To prove (B4), we observe that for families where e s = 1, because of the factorization 
( 3-71) , for any values of tt and e g , the locus B is the unit circle \q — 1| = 1. The result on 
chromatic zeros follows immediately from ( |3.7|) and the fact that the zeros of the chromatic 



polynomial of the circuit graph lie exactly on the circle \q — 1| = 1 |l5 . 

To prove property (B5), we first observe that since p is even, the left-hand side of eq. 
( p.lOj ) has the value 2 for q = 2, i.e., for a = 1. The condition that p is even means that either 



both ei and e 2 are even or both e\ and e 2 are odd. In both of these cases, the right-hand 
side of eq. ( |3.10|) is also equal to 2 for q = 2, so that q = 2 is a solution of this equation. 
Furthermore one easily checks that (|3.10|) has no solution for larger real q. This yields (B5). 



For property (B6) pertaining to odd p, we recall that the case e s = 1 is already covered 
by (B4); for e s > 2, we first note that since p is odd, one of (ei, e 2 ) is odd and the other is 



5 In the technical terminology of algebraic geometry [23], a multiple point on an algebraic curve 
is a point where several branches of the curve intersect. If all rij of the branches have different 
tangents, the multiple point is said to have index raj. 
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even; these are denoted e even and e dd- We next observe that for odd p = 2k + 1, D p contains 
a factor (q — 2) times a polynomial: 



where 



D2k+i — (q — 2)F 2 fc + i 



F 2 k+i = k at q = 2 



Now let us consider the limit q — > 2. From eqs. ( [3.16| ) and (|3.17|) , we have 

(q-2)(p-l) 



Furthermore, 



«2 



_(_l)e S+ l (g _ 2)(1+ee 



as q 



eodd) as g -> 2 



Hence for the chromatic polynomial 

'(g-2)(p 
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Since 



1 1 "I - ^even ^odd I 



p — 1 



1 + 



< 1 



(3.16) 



(3.17) 



(3.18) 



(3.19) 



(3.20) 



(3.21) 



for the relevant case e s > 2 considered here, it follows that, taking the m — > oo limit first, 
before the g — > 2 limit, the contribution of the second term in the square brackets, i.e., the 
ci2 term, in eq. (|3.20| ), relative to that of the first, i.e., a± term, goes to zero. Recalling our 
order of limits ( |1.2|) , it follows that at q = 2, the function W is determined completely by 
the ai term. Hence, although ai and a 2 both vanish at q = 2, this point is not on the locus 
B since W is determined only by a\ at this point. This shows that the condition of the 
degeneracy of terms |di(g)| = |a2(<z)| is necessary in order that the point q lie on the locus 
B, but, as noted in connection with eq. (|3.9| ), a further necessary (and sufficient) condition 
is that this equality in magnitudes holds where a\ and a 2 are both nonzero leading terms. 
Since W is determined by a± at q — 2, the same as in region R\ (see below), this shows that 
q c < 2, which was to be proved. 

The other property in (B6) follows in a similar manner from the solution to eq. ( 3.10Q . 
The fact that q c — > 2, i.e., a c — > 1, where a c = q c — 1, is clear from (|3.10|) since if a > 1 
(a < 1), the left-hand side diverges (vanishes) as e g — > oo, precluding equality with the 
right-hand side. 
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We observe several other interesting features of the locus B. For e g = and e s > 2 
and for certain values of (ei, €2)-, this locus, as an algebraic curve, may contain one or more 
multiple point(s), which are of index 2, i.e., they involve a crossing of four curves forming 
two branches of B. A particularly simple set of families for which this occurs is the set where 
ei = e<i = p/2 = ei2- By an analysis of eq. ( p.lOj ), we find that for this set there are 

AW = eia - 1 (3.22) 

multiple points, which occur at 

qj = 1 + e ^ (3.23) 

where 

% = ±— , J = 0,1,...,^-^ for e 12 even (3.24) 

Cl2 * 

a ,(2j + l)vr . ei 2 -3 

0j = ± , j = 0,1,..., — - — for ei2 odd (3.25) 

ei2 ^ 

For example, in Fig. (@)(a) for (ei, e 2 , e g ) = (2, 2, 0) one sees that there is one such multiple 
point, located at q = q c = 2, while in Fig. @(a) for (ei,e2,e g ) = (3,3,0) there are 2 
complex-conjugate multiple points, located at q = 1 + e ±Mr / 3 . Fig. (|6]) illustrates another 
case with a multiple point, namely (ei,e2,e g ) = (2,4,0), for which the multiple point is at 
q = q c = 2. In Figs. 2-3 and 6-9 one sees that the multiple points that are present for e g = 
disappear when e g > 1. Figs. 4-5 illustrate a family with p odd, viz., (61,62) = (2,3) for 
which B has no multiple points. 

In the cases where B contains multiple point (s) for e g = 0, which disappear when e g > 1, 
each such disappearance is accompanied by the appearance of a new disconnected component 
of B. As 6 g increases, the sizes of the regions bounded by these components of B decrease, 
and as e g — > 00, the regions shrink to points. For example, in Figs. (0) and (§) one sees that 
for the case (2, 2, e g ), as e g increases from e g = to e g > 1, the multiple point at q = q c = 2 
that was present for e g = disappears and there appears a disconnected component of B 
defining a new region, so that the number N comp , of (separate, disconnected) components of 
B increases from 1 to 2. As e g — > 00, this region shrinks to a point at g = 3/2. Similarly, 
in Figs. (H) and (§) one sees that for the case (ei,e 2 ) = (3,3), as e g increases from to 
e g > 1, the complex-conjugate multiple points that were present for e g = at q = 1 + e ±J7r / 3 
disappear and N comp , increases from 1 to 3. As e g — » 00, the inner regions shrink to points 
at q ~ 3/2 ± i/2. 
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Another feature that one observes in the results shown in the figures is that these loci 



B do not have any support for Re(q) < 0. This is in accord with our conjecture [19,22 



that global circuits are a necessary but not sufficient condition for B to have support for 
Re{q) < 0. 

For W([\im m _ 00 G eite2i e gjm \,q), which henceforth will be abbreviated simply as W, we 
have the following general results. In region R±, 

W = (ai)^^r = [( q - ly^Dp}^^ for q E R x (3.26) 

We note that the function on the right-hand side of eq. ( |3.26|) is the same as that for the 
open chain, (|2.10|) , although in the latter case, the result holds throughout the full q plane 
whereas for the cyclic chain, the result holds only in region R\. One can approach the origin, 
q = from the left staying within the region Ri, and in this limit, using the property ( |3.15| ), 
we have 

W = (p- l)f^=r for q -> , qeRi (3.27) 

Let us denote region R 2 as the one containing the point q = 1; this region is contiguous to 
Ri at the origin, q = 0. In region R 2 , a 2 is the dominant term, so that 

\W\ = \a 2 \^+^ 1 for qeR 2 (3.28) 

In particular, for families of graphs with e s = 1, where the factorization ( |3.7| ) holds, 

\W\ = {Dpl^^ 1 for e s = 1 and q G R 2 (3.29) 

For cases where there are additional regions Rj enclosed within R 2 , 

i 

| W| = |ai | p+^-i for q e Rj enclosed within R 2 (3.30) 

For e g = these additional regions are contiguous with Ri at multiple points on B. 

Concerning values of W at special points, we observe that since q c < 2, the expression 
( |3.2(j| ) always holds at q — 2, so that, as a consequence of the property (fTT§), 

W(q = 2) = for odd p (3.31) 

Note that this does not follow just from the property that the chromatic polynomial P 
vanishes at q = 2 for odd p, but requires the stronger property that P contains the factor 
(q — 2) u with v oc n as n — > oo. In contrast, for even p = 2k, using the property 

D 2k = 1 for q = 2 (3.32) 
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we have 



W(q = 2) = 1 for even p (3.33) 
As one approaches q = from within region R 2 , using ( 3.15Q , it follows that 



\W\ = \p- l\p+°9-* for g^O, qER 2 (3.34) 

in accord with ( |3.27|) and the equality of the magnitude of \W\ across any point on B. 

In Figs. (@)-@ we have also shown the chromatic zeros of G eit e 2)e m for large finite values 
of m to compare with the m — > oo accumulation set B. We observe first that the chromatic 
zeros lie quite close to the asymptotic locus B. Second, the density of chromatic zeros 
appears to vanish at multiple points, just as was the case for multiple points on complex- 
temperature phase boundaries for which we carried out exact calculations previously ||24|| . 
As one increases e g above 0, so that new components of B pinch off, one still sees a remnant 
of the reduction of density on the sides facing the location of the multiple point that had 
existed for e g = 0. 
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-2.00 1 1 1 1 1 1 1 1 1 1 1 1 

-1.00 0.00 1.00 2.00 3.00 4.00 

Re(q) 

FIG. 2. Boundary B in the q plane for W function for linv^oo G e i,e 2 ,eg,m with (ei, e 2 , e g ) — (a) (2,2,0), 
(b) (2,2,1). Chromatic zeros for m = 14 (i.e., n = 42 and n — 56 for (a) and (b)) are shown for comparison. 
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-2.00 1 1 1 1 1 1 1 1 1 1 1 1 

-1.00 0.00 1.00 2.00 3.00 4.00 

Re(q) 

FIG. 3. Boundary B in the q plane for W function for linim^oo G e i,e 2 ,eg,m with (e\,e2,e g ) = (a) (2,2,2), 
(b) (2,2,3). Chromatic zeros for m = 10 (i.e., n = 50 and n — 60 for (a) and (b)) are shown for comparison. 
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2.00 



1.00 



f 1 0.00 



-1.00 




(a) 



-2.00 ' ■ 1 1 1 1 1 1 1 ■ 1 — 

-2.00 -1.00 0.00 1.00 2.00 3.00 

Re(q) 

2.00 i ■ , ■ , ■ , ■ , ■ , — 



1.00 



•f 0.00 



-1.00 




(b) 



-2.00 ' 1 1 1 1 1 1 1 1 1 1 1 

-2.00 -1.00 0.00 1.00 2.00 3.00 

Re(q) 

FIG. 4. Boundary B in the q plane for W function for linv^oo G ei ,e2,e 9 ,m with (ei, e 2 , e g ) = (a) (2,3,0), 



(b) (2,3,1). Chromatic zeros for m = 14 (i.e., n = 56 and n = 70 for (a) and (b)) are shown for comparison. 
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-2.00 1 ■ 1 1 1 1 1 1 1 ■ 1 — 

-2.00 -1.00 0.00 1.00 2.00 3.00 

Re(q) 

2.00 i ■ , ■ , ■ , ■ , ■ , — 



(b) 



1.00 



•f 0.00 



-1.00 




-2.00 ' 1 1 1 1 1 1 1 1 1 1 1 

-2.00 -1.00 0.00 1.00 2.00 3.00 

Re(q) 

FIG. 5. Boundary B in the q plane for W function for linim^oo G ei>e2ie m with (e\,e2,e g ) = (a) (2,3,2), 
(b) (2,3,3). Chromatic zeros for m = 10 (i.e., n = 60 and n = 70 for (a) and (b)) are shown for comparison. 
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IV. CONCLUSIONS 



In this paper we have given exact expressions for the chromatic polynomial P(G, q) and 
the resultant exponent of the ground state entropy for the Potts antiferromagnet in the 
n — > oo limit, W({G}, q) for families of cyclic polygon chain graphs G = G eit e 2 ^ Stm . We have 
studied several types of limits yielding n — > oo, viz., , L £g , L p , and L m , the last of which 
yielded the most interesting results. From these calculations we have derived the continuous 
locus, £>, of nonanalyticities of W , which is also the accumulation set of the zeros of the 
chromatic polynomial in the n — > oo limit. Various properties of this locus were proved. A 
comparison with the results for the open polygon chain graph shows the important effect of 
global circuits. The results of the present study agree with and extend our earlier inferences 
concerning the locus £>, in particular, that a sufficient condition that in the n — > oo limit the 
locus B separates the q plane into two or more regions is that the graph has a global circuit 
with lim^oo £ g _ c _ = oo. 

This research was supported in part by the NSF grant PHY-97-22101. 



V. APPENDIX 

In this appendix we present results for the families of cyclic polygon chain graphs with 
(ei,e 2 ) = (2,4) and (3,3), and, for each case, < e g < 3. 
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-1.40 1 1 1 1 1 1 ' 

— 1 .00 0.00 1.00 2.00 

Re(q) 




-1.40 1 1 1 1 1 1 1 1 1 

-1.00 0.00 1.00 2.00 

Re(q) 

FIG. 6. Boundary B in the q plane for W function for linim^oo G ei! e 2i e g , m with (ei, e 2 , e g ) — (a) (2,4,0), 
(b) (2,4,1). Chromatic zeros for m = 10 (i.e., n = 50 and n — 60 for (a) and (b)) are shown for comparison. 
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FIG. 7. Boundary B in the q plane for W function for lim™-^ Gei,e 2 ,e 9 ,m with (ei, e 2 , e g ) = (a) (2,4,2), 
(b) (2,4,3). Chromatic zeros for (a) m = 8 (i.e., n = 56) and (b) m — 6 (i.e., n = 48) are shown for 
comparison. 
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2.0 



1.0 



■f 0.0 



-1.0 




-2.0 ' ■ 1 1 1 1 1 ■ 1 ■ ' 

— 2.0 -1.0 0.0 1.0 2.0 3.0 

Re(q) 

2.0 i ■ , ■ , ■ , ■ , ■ ^- 



1.0 



f 0.0 



-1.0 




(b) 



-2.0 ' 1 1 1 1 1 1 1 1 1 1 1 

-2.0 -1.0 0.0 1.0 2.0 3.0 

Re(q) 

FIG. 8. Boundary B in the q plane for W function for linv^oo G e i,e 2 ,eg,m with (e\,e2,e g ) = (a) (3,3,0), 
(b) (3,3,1). Chromatic zeros for m = 10 (i.e., n = 50 and n — 60 for (a) and (b)) are shown for comparison. 



23 



2.0 



(a) 




-2.0 -1.0 0.0 1.0 2.0 3.0 

Re(q) 

FIG. 9. Boundary B in the q plane for W function for linim^oo G ei , e2 ,e 9 ,m with (e\,e2,e g ) — (a) (3,3,2), 
(b) (3,3,3). Chromatic zeros for m = 8 (i.e., n = 56) and n = 64 for (a) and (b)) are shown for comparison. 
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